All inseparable two-mode Gaussian continuous variable states are distillable 
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We show that all entangled Gaussian states of two infinite dimensional systems can be distilled to 
maximally entangled states in finite dimensions. The distillation protocol involves local squeezing 
operations, local homodyne measurements with ancilla systems prepared in coherent states, and 
local joint measurements of the total number operator of two copies of the state. 
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The existence of pure entangled states of two or more 
systems entails the possibility of finding new applications 
of Quantum Mechanics, in particular in the fields of com- 
putation and communication. In the latter case, sev- 
eral protocols were devised that use entangled pairs of 
qubits (quantum systems described in terms of a Hilbert 
space of dimension two) to achieve secure secret com- 
munication Q| . These protocols have been successfully 
implemented experimentally by using pairs of photons 
generated via parametric down conversion H. More re- 
cently, these protocols have been generalized to the case 
of infinite dimensional Hilbert spaces, i.e., the so called 
continuous variable (CV) systems JH . The successful ex- 
perimental realization of teleportation || using light 
beams in squeezed states indicates that those (or more 
efficient) protocols may be very soon implemented in the 
laboratory. It is clear that with light beams in squeezed 
states (containing many photons) one could increase the 
transmission rate in quantum communication compared 
to the systems realized today. 

In practice, however, systems are exposed to interac- 
tions with the environment. These interactions trans- 
form pure into mixed states, which may not be useful for 
quantum communication. Fortunately, there exist meth- 
ods to recover pure entangled states from mixed ones in 
certain situations. These processes are called entangle- 
ment distillation (or purification) [tJ , and consist of lo- 
cal operations and classical communication transforming 
several copies of a mixed entangled state into one (ap- 
proximately) pure entangled state which is then useful 
for quantum communication. In fact, using this method 
in the appropriate way one can construct quantum re- 
peaters H that should allow quantum communication 
over arbitrarily long distances. Unfortunately, it is not 
known, in general, which mixed states p are distillable. 
At the moment we only have necessary or sufficient cri- 
teria for distillability, but not both. First, obviously the 
state p must be inseparable (entangled) for it to be distil- 
lable. Moreover, as shown by Horodecki et al. Q, there 
exists a stronger necessary condition, namely that p must 
have non-positive partial transpose. In fact, surprisingly 
enough there are entangled states which are not distil- 
lable since their density matrices remain positive under 
partial transposition JToj ] . Furthermore, there is evidence 
that this condition is not sufficient since there seem to 
exist states that have non-positive partial transpose but 



that nevertheless are not distillable |ll|]. On the other 
hand, a useful sufficient criterion, the so-called reduction 
criterion, has been reported Jl^] . It states that if there 
exists some vector 1-0) such that 



\tv B p® 1 -p\ip) < 0. 



(1) 



then the state p is distillable. Here, tr# stands for the 
partial trace with respect to the second subsystem. An 
important aspect of this criterion is that if one can find a 
state 1-0} satisfying (|J), then one can explicitly construct 
a protocol to distill p. 

Although all these results can be easily generalized to 
CV systems, in that case it is very complicated (if not im- 
possible, in general) even to determine whether a mixed 
state is entangled or not. However, for a class of states 
comprising nearly all the states that can currently be gen- 
erated in the laboratory, a necessary and sufficient condi- 
tion for inseparability has been derived [jl3| Q (in fact, in 
JTij it is shown that the partial transposition condition is 
also sufficient). Those states are called Gaussian states, 
since their density matrices can be written as Gaussian 
functions of position- and momentum-like operators act- 
ing on the two infinite-dimensional Hilbert spaces of the 
two subsystems. As mentioned above, this necessary and 
sufficient condition for entanglement is, in general, not 
sufficient for a mixed state to be distillable. In this Let- 
ter, however, we show that for Gaussian states of two 
modes inseparability implies distillability: all entangled 
Gaussian states of two modes are shown to be distillable 
and thus useful for quantum communication. Moreover, 
our proof, which is based in part on the reduction cri- 
terion, provides an explicit protocol that accomplishes 
distillation for all those states. 

Let us consider a Gaussian state p of two systems A 
and B acting on a Hilbert space L 2 (R) ® L 2 (R), e.g. two 
modes of the electromagnetic field. (We will in the fol- 
lowing often use quantum optical terms, like "modes" or 
"beam splitters", because quantum optics currently of- 
fers the most promising setting for the realization of CV 
systems. The results of this paper, however, are valid for 
all CV systems.) For such systems, it is convenient to 
describe the state p by its characteristic function 111] 



X(£)=tr[p£>(£)]. 
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Here £ = (q a ,p a , q b ,p b ) G R 4 is a real vector and 

— e - i Y. k = a,b(lkX k +p k P k ) ^ 

where X a ^ b and P ajb are operators of system A and B, 
respectively, satisfying canonical commutation relations 
(Ti = 1). The characteristic function contains all the 
information about the state of the system (e.g., JL5[), 
that is, one can find p starting from \- Gaussian states 
are exactly those for which \ is a Gaussian function of £ 



x(0 



(4) 



where the correlation matrix M is a real symmetric ma- 
trix, the displacement d <E R 4 a real vector, and where 
(■,•} denotes the standard scalar product in R 4 . Thus, 
a Gaussian state is fully characterized by the correlation 
matrix M and the displacement vector d. 

Any Gaussian state of two modes can be transformed 
into what we called the standard form, using local unitary 
operations only 13, H. The corresponding characteris- 



tic function has displacement d = and the correlation 
matrix M has the simple form 



M = 



M A Mab 
M T AB M B 



(5) 



where 



M A = 



n a 
n a 



,Mb = 



n b 
n b 



Mab = 



K, 



(6) 



As mentioned in |13|] the local unitaries needed to 
achieve standard form are linear Bogoliubov transforma- 
tions, i.e., generated by Hamiltonians that are at most 
quadratic in the operators X a ^,P a ,b- For optical fields 
this means that they can be performed with beam split- 
ters, phase shifters and squeezers, i.e., currently available 
technology. 

The four real parameters (n a ,n b ,k x ,k p ) fully charac- 
terize a Gaussian state up to local linear Bogoliubov 
transformations (LLBT). They can be easily calculated 
from the four LLBT-invariants \M A \ ,\M B \, \M A b\, and 
\M\ as follows: 



\M A 



\Mb \ ,k x k p 



|M 



AB 



(n a n b - k 2 x )(n a n b - k 2 p ) = \M\, 



(7a) 



(7b) 



where \M\ denotes the determinant of M; without loss 
of generality we choose k x > \k p \. 

Note that not all values of these parameters are al- 
lowed, since the correlation matrix must correspond to 
the density operator of a physical state. It has been 
shown p7| that a given correlation matrix M corresponds 



to a physical state if and only if (iff) M > and —JMJ > 
M , where J is a 4 x 4 matrix called complex struc- 
ture, defined by J(q a ,p a ,qb,Pb) T = (-p a ,q a ,~Pb,qb) T ■ 
In terms of the parameters (^) we can reexpress these 
conditions as 

(n a n b - kl)(n a n b - k p ) + 1 > n\ + n\ + 2k x k p ,(8a) 
n a n b — k x > 1. (8b) 

On the other hand, starting from results of [Q it is 
easy to show that a Gaussian state is entangled iff the 
corresponding parameters satisfy 



(n a n b - k 2 x ){n a n b - k p ) + 1 < n\ + > 



2k x kp. (9) 



In the following we will prove that a Gaussian state 
is distillable iff its parameters fulfill (^|); that is, iff it 
is entangled. We will proceed as follows: first, we recall 
the proof of the reduction criterion in |{l2| and extend the 
result to infinite dimensions. The extended criterion pro- 
vides a sufficient condition for distillability of CV states. 
Secondly, we show for symmetric Gaussian states, i.e. 
states for which n a = n b = n, that the inseparability 
criterion (^|) is equivalent to the reduction criterion. So, 
all the symmetric inseparable states are distillable. Fi- 
nally, we show that all the states which are not symmetric 
(n a 7^ nb) can be brought into a symmetric form by using 
local operations and without disentangling them. 

We start by reviewing the proof of the reduction cri- 
terion for distillability and generalizing it to infinite di- 
mensions. Let a density matrix p and the pure state 
\ip) = J2 n m a nm \n) ® \m) fulfill the criterion (|j), where 
the vectors \n) form an orthonormal basis of TL. The 
coefficients a nm define a matrix A — [a nm ) satisfying 
AA^ = trs(|-0) Distillation of p is divided into 

three steps. The first is a filtering operation: The oper- 
ator AA* <8> 1 can be viewed as an element of a positive- 
operator- valued measure (POVM), which defines a gen- 
eralized measurement Jl§[ ]. Conditional on the measure- 
ment outcome corresponding to AA^ 1 we obtain the 
state 



(10) 



p = A* ® IpA ® l/tr(pAA* O 1), 

which still satisfies (Q) but now with \ip) = |<i>+) := 
^=5ZfcLil^i^)> the symmetric maximally entangled 
state of two TV-level systems. In this case, (Q) implies 
tr(p|$+) ($+|) > l/N. A state satisfying this inequal- 
ity can be distilled by a generalization of the protocol 
of Ref. |J, which consists of two steps: depolariza- 
tion and joint measurements. Applying an operation of 
the form U <8> U* (U unitary, randomly chosen) depo- 
larizes p, i.e. transforms it into a mixture of the max- 
imally entangled state (which is invariant under 
transformations of the form U ®U*) and the completely 
mixed state j^s-l; the overlap of p with 3>+) remains 



3 



unchanged. Taking two entangled pairs in this depolar- 
ized form, both A and B perform the generalized XOR 
gate XORjv : \k) \l) >— ► \k) \(l + k)modN) on their respec- 
tive systems. Then they both measure the state of their 
second system in the basis \k). The first pair is kept, 
iff they get the same result. The resulting state has a 
density matrix p' , which has a larger overlap with the 
maximally entangled state |$+) than the original p. It- 
erating the last two steps sufficiently often, the overlap 
between the resulting state and |$+) approaches 1, that 
is, we obtain a nearly maximally entangled state. 

The reduction criterion was proved in |l2j for finite 
dimensional systems. Since we want to apply it to CV 
states we have to show that it extends to infinite dimen- 
sions, i.e., that Ineq. (|j) implies distillability even for 
dimTi = co. Let {\k) : k = 0,1,...} be an orthonormal 
basis of H, let TL n — span{|0) , |1) , \n)}, and let p be a 
density matrix on 7i ®H. Let £ (p) = trap® 1 — p be the 
map occuring on the lefthand side of |l]). Assume that 
3 € H, e > such that {ip\ £(p) < -e < 0. Since 
Pn = Pn n pPn n converges weakly to p, there is N > 
such that (ip\ S(p n ) < -e/2 for all n > N. Thus p 
can be projected by local operations to a distillable state 
Pn and is therefore itself distillable. ■ 

Now we specialize to Gaussian states. If both the states 
p and ip occuring in (|l|) are Gaussian with displacements 
dp = d^p = and correlation matrices M p , respec- 
tively, then ([!]) takes the form 

2 [|M trBP + M trf ^|]- 1/2 - 4 [\M„ + M it \r 1/2 < 0. (11) 

Let p be a symmetric Gaussian state in standard 
form and %jj the pure two- mode squeezed state \tp) — 
— tt- Y)„ tanh" r Inn). This is also a Gaussian state and 

cosh r £-~* n I J— 

the four parameters (|7|) are n a = rib = coshr, k x — —k p = 
sinhr. In the limit of large r (keeping only the leading 
terms in e r ) this becomes after some simple algebra 



(n + k x )(n — k p ) > 1. 



(12) 



For symmetric states, however, Ineq. (12) is implied 
by the inseparability criterion: for symmetric states, @ 



simplifies to 



k x k p — 1| < n(k x — k p ), which implies 



k x kp > 1 or < 1, proving that all 



Ineq. (@) both if n 
symmetric inseparable Gaussian states are distillable. 

If the state is not symmetric, it means that the reduced 
state at one of the two sides has larger entropy than the 
other. This suggests to let a pure state interact with 
the "hotter" side to cool it down. To do this without 
destroying the entanglement of p, we proceed as follows: 
p is transformed such that the correlation matrix of its 
Wigner function takes on its standard form with param- 
eters (N a , Nb, K x , Kp). The Wigner function is the sym- 
plectic Fourier transform of the characteristic function 
([|) and therefore also Gaussian for Gaussian states. The 
Wigner correlation matrix Mw is related to the (char- 
acteristic) correlation matrix by Mw — —JM~ X J [fill. 



One can formulate the conditions (^,^|) similarly in terms 
of the parameters (N a , Nb, K. x , K p ). Just as the form 
(§), the Wi gner standard form can be obtained by local 
squeezing operations. 

Now assume that Nb < N a , i.e., B is the hotter side. 
Take an ancilla mode in the vacuum state and couple it 
to B's mode by a beam splitter with transmittivity cos 2 9. 
After a homodyne measurement of the ancilla results a 
state p with Wigner correlation matrix M with 



- If c 2 N a + s 2 D x 
MA ~u\ c 2 N a + s 2 N a N b 



w _ 1 f N b 

MB ~v\U (c 2 N b + s 2 )v 



Mai 



cK x 




cKpV 



where the abbreviations c = cos#, s = sin#, v = s 2 Nb + 
c 2 , and D x p — N a Nb — K x were used. The condition 
for symmetry, \Ma\ — \Mb\, requires 



tan 2 6 



N b - D x N a 



(13) 



Checking @ for M one easily sees that the sign of the 
lefthand side does not change; therefore the transformed 
state is inseparable iff the original one was inseparable. It 
remains to be shown that there always exists a 6 to satisfy 
(|l3|), i.e., that the right hand side of Eq. ( p^ ) is positive. 
The numerator is positive since we have chosen Nb < N a , 
the denominator is positive for all states since the second 
part of condition (^) implies that (N a — D x Nb) > and 
the first part assures that (N a — D x Nb)(Nb — D p N a ) > 
(N a K x + NbKp) 2 > 0, hence all Gaussian states can be 
symmetrized this way. But since every Gaussian state 
can be brought into Wigner standard form by local uni- 
taries, this completes the proof. ■ 

Moreover, this implies that the protocol of Ref. 
can in principle be used to obtain maximally entangled 
states in a finite dimensional Hilbert space from any given 
inseparable Gaussian two-mode state. However, a prac- 
tical distillation protocol for all Gaussian states remains 
to be found. The realistic proposals for Gaussian states 
have only been shown to work for pure states p(| or in 
the limit of small or one-sided noise H] . It is worth not- 
ing, though, that major parts of the universal protocol 
of p2| are in fact implemented by the procedure of Duan 
et al. pl| . Let us briefly discuss the steps for distillation 
of an inseparable Gaussian state p: 

(i) remove displacements d ^ 0; symmetrize p if nec- 
essary; bring the symmetric state into standard form. 
These steps require the local use of beam splitters, one- 
mode squeezers, ancilla systems in coherent states, and 



4 



a homodyne measurement. 

(ii) for a state in symmetric standard form the filter- 
ing operation ( JlOj ) is unnecessary, since then p already 
satisfies Ineq. (J^) with state oc lim^i J2k ^° 1^) 1^) 
(in the photon number basis). This gives a nrn — 
lim A ^i \ n+m 5 nm , hence A = 1. 

(iii) depolarize the state; in |p| this step is not neces- 
sary; we do not know how to completely depolarize an 
arbitrary state quantum optically. However, since all this 
operation does is to increase our ignorance of the system 
it is not necessary for distillation but could for any given 
p be replaced by some specific unitary that would work 
just as well or even better. 

(iv) This is the most important step of the distillation 
protocol. A bilocal XOR is used to mutually entangle 
two entangled pairs. A subsequent measurement selects 
a distilled subenscmblc. This operation is implemented 
in |2l[ by the total photon number measurement at both 
sides: The state conditional on both A and B obtaining 
the same result N differs only by a local unitary trans- 
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forn iation (namely \ n,N — nf- 



obt; 



N)) from the one 



incd by fir3t projecting bilocally to the N I 1 di - 



mensional subspace Hn+i {p — ► Pn+i), performing the 
bilocal XORat + i and finally measuring the target sys- 
tem. As shown before, for a sufficiently large value of N, 
the truncated state pn+i is distillable and then step (iv) 
leads to a state with larger overlap with 
Iteration of the steps (i)-(iv) leads with finite probability 
arbitrarily close to a maximally entangled state in some 
finite dimensional space. 

In conclusion we have shown that all inseparable Gaus- 
sian states of two modes are distillable into maximally en- 
tangled states in some finite dimensional Hilbert space. 
We have described a protocol that accomplishes this and 
discussed its relation to the quantum optical distillation 
procedure of Ref. piJl , showing that major parts of the 
protocol can be implemented by simple quantum opti- 
cal means. A number of open questions remain: firstly, 
it is not clear how to implement the depolarizing oper- 
ation quantum optically. It will be further investigated 
whether this step is necessary at all for Gaussian states 
and how it might be performed quantum optically. To 
iterate the protocol, the same question has to be an- 
swered for the states obtained by (repeatedly) perform- 
ing steps (i)- (iv) on a Gaussian state. Note also that 
the protocol [£1[ works better, when the total photon 
number measurement is performed on more than two en- 
tangled pairs simultaneously (indeed being optimal for 
pure states asymptotically). If using many modes simul- 
taneously improves the protocol for mixed states, too, 
remains to be seen. Furthermore, it would be interesting 
to see whether the equivalence of inseparability and dis- 
tillability also holds for Gaussian states of tua + Wb > 2 
modes or whether one may find bound entangled Gaus- 
sian states of bipartite systems of three or more modes. 
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